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ABSTRACT 

The fast variability of energetic TeV photons from the center of M87 has been detected, offering a new clue to 
estimate spins of supermassive black holes (SMBHs). We extend the study of Wang et al. (2008) by including 
all of the general relativistic effects. We numerically solve the full set of relativistic hydrodynamical equations 
of the radiatively inefficient accretion flows (RIAFs) and then obtain the radiation fields around the black hole. 
The optical depth of the radiation fields to TeV photons due to pair productions is calculated in the Kerr metric. 
We find that the optical depth strongly depends on (1) accretion rates as t TcV oc M 2 ,5 ~ 5 °; (2) black hole spins; 
and (3) location of the TeV source. Jointly considering the optical depth and the spectral energy distribution 
radiated from the RIAFs, the strong degeneration of the spin with the other free parameters in the RIAF model 
can be largely relaxed. We apply the present model to M87, wherein the RIAFs are expected to be at work, 
and find that the minimum specific angular momentum of the hole is a ~ 0.8. The present methodology is 
applicable to M87-like sources with future detection of TeV emissions to constrain the spins of SMBHs. 
Subject headings: black hole physics — galaxies: individual (M87) 



1. INTRODUCTION 

Astrophysical black holes can be simply depicted by two 
parameters: mass and angular momentum. Masses of super- 
massive black holes (SMBHs) are relatively easier to estimate 
since their gravitational influences spread over the large-scale 
space approached by Newtonian mechanics. SMBH masses 
in active galactic nuclei or galactic centers can be measured 
by several different methods through stellar or gas dynam- 
ics, or reverberation mapping of emission lines (Kaspi et al. 
2000; Kormendy & Gebhardt 2001; Ho 2008). However, 
spins are more elusive to estimate because their general rela- 
tivistic (GR) effects only appear significantly in the innermost 
region close to the black hole, typically within ^20 gravita- 
tional radii and eventually disappear outward then. Spatially 
resolving the region around ~10 or a few gravitational radii 
is in principle plausible through radio observations for the 
Galactic center (Shen et al. 2005; Doeleman et al. 2008), 
but it is still a limit challenge to current radio astronomy for 
estimation of the spins of the extragalactic SMBHs. 

Hitherto, only a few Seyfert galaxies show relativistically 
broadened, highly redshifted iron Ka emission lines, which 
are interpreted most plausibly by the collective effects of the 
Doppler motion of the fluids from which the intrinsic nar- 
row emission lines originate, the gravitational redshift, and 
the gravitational lensing around a rapidly rotating black hole 
(see, e.g., Miller 2007 and references therein). The most con- 
vincing evidence for effects of spins is found in the MCG-6- 
30-15 from XMM-Newton observations (Fabian et al. 2002). 
On the other hand, cosmic X-ray background radiation sug- 
gests that most SMBHs are fast rotating (Elvis et al. 2002) 
and the Soltan's argument applied to large samples of quasars 
and galaxies indicates a similar conclusion (Fabian & Iwa- 
sawa 1999; Yu & Tremaine 2002; Wang et al. 2006). 

TeV photons suffering from attenuation by pair productions 
are expected to explore the radiation fields in the vicinity of 
SMBHs and thus constrain their spins (Wang et al. 2008, 



hereafter Paper I). Rapid variability of TeV emission (at a 
timescale of 2 days) has been discovered by the HESS (High 
Energy Stereoscopic System) collaboration in the famous ra- 
dio galaxy M87 (Aharonian et al. 2006). Interestingly, the 
TeV emission does not originate from a relativistic jet (Aha- 
ronian et al. 2006), in contrast to cases of blazars (Blandford 
& Levinson 1995; Levinson 2006). A novel mechanism for 
TeV emission around the horizon of a spinning black hole by 
a magnetospheric pulsar-like process was originally proposed 
by Levinson (2000), and had been applied to M87 (Neronov 
& Aharonian 2007, but see Rieger & Aharonian 2008). In 
such a context, the TeV photons are able to serve as a probe 
of spins which determine the density of the radiation fields. 
Wang et al. (2008) found that the SMBH in M87 should have 
specific angular momentum a > 0.65 in order to allow the TeV 
photons to escape from the innermost region of the radiation 
fields from the accretion disk. However, their calculations are 
based on the self-similar solutions of ADAFs in Newtonian 
approximation and the GR effects on the disk structures and 
the propagation of photons are not included. 

The main goal of this paper is to extend the study of Paper I 
by including all the GR effects. In Section 2, we introduce the 
GR RIAF model. Section 3 gives the description of relativistic 
optics and the detailed formulae used for calculation of optical 
depth to TeV photons. The results are presented in Section 4 
with exhaustive investigation of parameter dependence of the 
optical depth. We then apply the results to M87 in Section 
5. Discussions and a summary are given in Section 6 and 7, 
respectively. Mathematical preliminaries are presented in the 
Appendix. 

2. GENERAL RELATIVISTIC RIAF MODEL 

We follow the work of Manmoto (2000) to construct 
the fully relativistic hydrodynamical equations. GR nota- 
tions are given in the Appendix. Cylindrical coordinates 
(t,R,cf),z) are used to describe accretion flows by expanding 
the Boyer-Lindquist coordinates around the equatorial plane 
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up to (z/R)° terms. All the basic equations can be derived 
from the conservation laws of mass, momentum, and energy 
under the four- dimensional space. For a clarification, we list 
the necessary equations as follows. 
The continuity equation reads 



M = -27rA 1/2 E 7 r V, 



(1) 



where V is the radial velocity of the flows with respect to the 
corotating reference frame (CRF, see the Appendix) with j r 
as its Lorentz factor (e.g., Abramowicz et al. 1996; Gammie 
& Popham 1998; Yuan et al. 2009), M is the mass accretion 
rate, and Eo is the surface density of accretion flows. The 
conservation of momentum gives two equations. The radial 
component is 

1 dW ll^M. (Q-Q+xn-ft-) 



dV 



7 Z V = — 

lr dR M S dR 



R 4 A 



(2) 



where 7^ is the Lorentz factor of the azimuthal velocity of the 
CRF with respect to the locally nonrotating frame (LNRF), 
the height-integrated total pressure W consists of gas pressure 
and magnetic pressure, and the relativistic enthalpy [i is writ- 
ten as 



fl=l + : 



W 



a, + - 



1 



+ a e + 



W 



Here we express the functions a, and a e as 



1 



7,-1 
1 



a e = 



2(1 ~/3 d ) 
ft ' 
2(1 -/3 d ) 



with the adiabatic indices of gas 7; and j e as 

"3/r 3 (iM)+tfi(iM) 



7; = l+0; 



4K 2 (l/0d 
3K 3 (l/6 e )+K 1 (l/6 e ) 



-1 



(3) 



(4) 



(5) 



(6) 



(7) 



4K 2 (i/e e ) 

where the global parameter /3d is the ratio of the gas pres- 
sure to the total pressure, = kTj/m p c 2 and 9 e = kT e /m e c 2 are 
the dimensionless temperatures of ions (7}) and electrons (T e ), 
respectively, k is the Boltzmann constant, m p and m e are pro- 
ton's and electron's mass, and K\,Ki, and ^3 are the modified 
Bessel functions. The azimuthal component of the conserva- 
tion of momentum is 



M(t-lin) = 2TTRW$, 



(8) 



where £ is the specific angular momentum of accretion flows, 
£i n is the swallowed by the central black hole, and 

A 3/2 A l/2 7 3 

-W, (9) 



W« = a d - 



R 6 



where ay is the viscosity parameter. The conservation of en- 
ergy for ions and electrons gives 



. ds 

-MT i -±=-2na i {\-5)W 
dR 

dR 



iIa 2 dn 



R 6 
dO 



dR 



■27r7?A„ 



(10) 



R* ar^^-F-), (11) 

where 5, (s e ) is the specific entropy of ions (electrons), A ie 
is the energy transfer rate from the ions to the electrons per 



unit surface area, and F~ is the radiative cooling per unit sur- 
face area. <5 is the fraction of the viscous dissipation which 
heats the electrons. To complete the set of equations, thermal 
dynamical relations among the entropy, the height-integrated 
pressure and the surface density for ions and electrons are 
given by 



Tidsi = 



T e ds e = 



W t 1 

w„ 1 



[d In W; - T, d m £ + (r, - 1 )d In fl] , (12) 



[dlnW e -r e dlnS + (r e -l)dln/?], (13) 



s e rw 

where the effective adiabatic indices are 

din a, 



r, = i + 



1 + 



1 + - 



dln7] 



1 -1 



1 + - 



dlna e 
dlnT, 



(14) 



(15) 



The vertical scale height H of the accretion flows is taken as 
(see also Abramowicz et al. 1997) 



So £ 2 -a 2 (e 2 -Li 2 y 



(16) 



where e is the specific energy. The simplification 
dln///dln/? = 1 (Manmoto 2000) is used in deriving the en- 
ergy equations. 

We take into account three processes of the radiative cool- 
ing, i.e., the synchrotron radiation, the bremsstrahlung, and 
the multi-Comptonization of soft photons. The general de- 
scription of cooling processes and relevant formulae have 
been presented by Narayan & Yi (1995) and Manmoto (2000) 
in a more handy way. In our calculations we completely make 
use of the program of the Comptonization given by Coppi & 
Blandford(1990). 

The outer boundary conditions of the GR RIAF model are 
imposed at a radius r out = R ou t/R g = 2 x 10 4 : 



fl = 0m K ; Ti = T e = 0.1T yi 



(17) 



where R g = GM./c 2 and r vir is the virial temperature defined 

as 



7W = (7-1) 



GM.m p 

kR 



1.1(7-1) x 10V K 



(18) 



with the adiabatic index 7 = 5/3 and = r/10 4 . Electrons 
and ions have the same temperature at the outer boundary. As 
shown by Manmoto et al. (1997), the outer boundary condi- 
tions have little influence on the structures of the GR RIAFs. 
We confirmed this effect and thus fix the outer boundary con- 
ditions as Equation (17). 

The choice of the eigenvalue l- m should satisfy the condi- 
tion that the flows pass through the sonic point smoothly. The 
inner edge of accretion flows is elusive since it depends on the 
accretion rate, and is nonaxisymmetric and time-variable (see, 
e.g., Krolik & Hawley 2002). We use the horizon of the black 
hole as the inner edge for a conserved influence of the spins 
on the optical depth of TeV photons. We show this effect of 
inner edge on uncertainties of spins in the following sections 
in details. 
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3. GENERAL RELATIVISTIC OPTICS 

3.1. Ray tracing method 

We assume that the accretion flows are axisymmetric, so 
the radiation fields have the same symmetry, and only the mo- 
tions on the (R-0) plane are required in our calculations. 
The general trajectories of photons in Kerr metric are de- 
scribed by three constants of motion (Bardeen et al. 1972). In 
terms of the photon's four-momentum, the conserved quan- 
tities are E = —p t , the total energy at infinity; L = p^, the 
component of angular momentum parallel to symmetry axis 
and q = Pq + cos 2 9a 2 p 2 + p\ cot 2 9. By defining A = p^/E and 

Q = q/E 2 , the four- momentum of photons can be rewritten as 



where 

K(r) = r 4 + (a 2 - 



A 2 -Q)r 2 + 2[Q + (A- 



-X z cot £ e 



2/ 



a 2 Q, 



E, 
(19) 

(20) 
(21) 



0(9) = Q + a 2 cos 2 9- 

where r = R/R g . 

Basically, the trajectory of a photon on the (R-9) plane is 
governed by the geodesic equations 

,7 = ±[^ = ±f^L= (22) 

where & is the affine parameter and the ± signs represent the 
increment (+) or decrement (-) of r and 9 coordinates along 
the trajectory. The reader is recommended to refer to Rauch 
& Blandford (1994), Cadez et al. (1998), and Li et al. (2005) 
for a comprehensive description of the analytic solutions of 
Equation (22). Given two constants of motion, A and Q, and 
the spin parameter a, the trajectory of a photon is uniquely 
determined. We neglect the photons with trajectories return- 
ing to the accretion disk but may encounter the TeV photons 
along their path because of rare probability of occurrence. 

3.2. The observed emergent spectrum 

The specific flux density observed by a remote observer is 
expressed as 

F u . = f I Vo dft , (23) 



where I Vo is the specific intensity as a function of frequency 
v in the observer's frame and df2 D is the element of the solid 
angle subtended by the image of the accretion disk on the ob- 
server's sky. In a common way, the image of the disk can be 
described by two impact parameters a and j3 (see Figure 1), 
which respectively represent the displacement of the image 
perpendicular to the projection of the rotation of the black 
hole on the sky and the displacement parallel to the projection 
of the axis (Li et al. 2005). Applying the invariant l v jir along 
the path of a photon (Rybicki & Lightman 1979, p. 146), we 
have 

dad/3 (24) 



F,, = 



D 2 



where v e is the frequency of the photons in the local rest frame 
(LRF) of the accretion flows, g = v /v e is the redshift factor 
of the photons (see Section 3.3), 7„ e is the specific intensity 
of the disk radiation at radius r e , and D is the distance of the 
black hole from the observer. Then the observed emergent 
spectrum is 



L„ =4irD l F„ o =4tt / g-%£a&p. 



h 



(25) 



For each (a,/3) set, the constants of motion A and Q can be 
expressed as (Li et al. 2005) 



A = -asin0 



Q = /3 2 + (a 2 - fl 2 )cos 2 e obs , 



(26) 



where the viewing angle o b s is the inclination between the 
rotation axis of the accretion disk and the direction to the ob- 
server. Using the ray tracing method we can locate the emis- 
sion place r e in the accretion disk for the photons which reach 
the observer's sky at point (a, (3), and, therefore, obtain the ra- 
diation intensity I Vt from the GR RIAF model. The observed 
emergent spectrum is obtained by integrating I Ve over the ob- 
server's sky as shown in Equation (25). 

3.3. Optical depth to TeV photons 

Figure 1 shows the geometric scheme for calculations of 
the optical depth to TeV photons which emanate from the po- 
sition P(/? TeV ,6Tev). The TeV photons unavoidably interact 
with the soft photons from the accretion disk when they are 
moving out [e.g., at the interacting point Q(r C7 9 c )], setting up 
strong constraints on the radiation fields of soft photons. We 
delineate the interactions in the LNRF (see Figure 2). Firstly 
we divide the solid angle at the interacting point in the LNRF 
into numerous elements (A9 s ,A(j) s ). For each solid angle el- 
ement at which the soft photons from the disk arrive, we can 
determine their constants of motion A s and Q s (see Equation 
(36)). Secondly we use the ray tracing method to trace the 
soft photons back to the accretion disk. We hence obtain the 
number density of soft photons. Lastly the optical depth is 
calculated by integrating all the soft photons from different 
directions along the trajectory of the TeV photons. 

Generally speaking, given the emission location of TeV 
photons (Rjev, OTev) and the viewing angle of an observer 
(8 bs), there is a large number of trajectories to the observer 
at infinity because of the arbitrariness of the (^-component. 
Since TeV photons along the shortest path intuitively suffer 
the least interactions and thus undergo the minimum optical 
depth, we only consider the optical depth for the shortest path 
in our calculations. We select the trajectory with a = as the 
shortest path, which represents the case that the TeV photons 
lie at the projection of the rotation axis on the image sky. Keep 
in mind that genuine emission position always has nonzero </>- 
component, i.e., «3>Tev 7^0, due to the drag of inertia caused by 
the rotation of black hole. 

In the LNRF, given the frequency of TeV photons v^ v at 
infinity, we can determine their corresponding frequency at 
the interacting point Q(r c ,9 c ) by redshift factor defined as 



8 TeV 



,TeV 



,TeV 



-if) 1 At 



pTeV 



£l/2 A l/2 X 



A 1 / 2 1-coX 



TeV 



(27) 



where the subscript r c or oo means their values are calculated 
at position (r c ,9 c ) or infinity and similarly hereinafter. In the 
same way, at the interacting point the frequency of soft pho- 
tons is given by the redshift factor in comparison with their 
frequency at the emanating place in the accretion disk, 



^(LNRRP^ 



e-Xl-wA,)!, 



< } (LRF)P U 



l-ftA s T 



1/2 



TJT 



(28) 

The directional angles of the TeV and the soft photons can 
be obtained from the projection of their four-momentum onto 
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observer's sky 



FIG. 1 . — Geometric scheme. The point P(7?tcV, ©Tcv) is the location of the TeV source and Q(r c , d c ) is the location of the interaction of TeV photons with the 
soft photons from the RIAFs. 




A* pTeV 



„M pTeV 
B (0 r A< 



A Te y SA'/ 2 
sin0 c A(l - wAxev) ' 



(31) 

and those of the soft photons are given by replacing the super- 
script (subscript) "TeV" with "s" in the above equations: 



e M P s 



FIG. 2. — Locally nonrotating frame at the interacting point Q(r c ,8 c ). 



the spatial directions of the LNRF, respectively. The direction 
cosines (a T(;V ,/3 Tt , v ,7 T( . v ) of the TeV photons in the LNRF are 
given by 



7 S = cos t 



j3, = sin 6 S cos < 



a s = sinf s sin< 



- P s 



p^ P* 



e <tre 



p^ P s 



- e' 1 P s 
c (0 P 



(32) 



(33) 



(34) 



Note that a 2 + /3 2 +7 2 = 1, just giving a s and f3 s , we can obtain 
A s and Q s from Equations (33) and (34). For simplicity, we 
define two denotations 



7tcV — CO *> ^TcV 



/^TeV — Sin $ Tc y COS y'i^ v 



M pTeV 



J 1 pTeV 



M pTeV 



±K{r) x l 2 
A l IH\-uu\ leN y 



8?-- 



a, sin 6 C A 



/3 s AV 2 (l-a;A s ) 



(29) 



EAV2 ' A [ /2 
then we can express A s and Q s by a s and /3 S as 



(35) 



A s = 



, Q % = m 2 -(acos6» c ) 2 + (AsCot6» c ) 2 . (36) 



pTeV 



A 1 /2(l- w A T ev)' 

(30) 



Having obtained A s and Q s , the ray tracing method is used to 
determine the emission location rj of the soft photons in the 
accretion disk. 
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The number density of the soft photons at the interacting 
point (r c ,0 c ) is 



CHI's cnv s 



(37) 



where h is the Planck's constant. Here we have applied the 
invariant I^/v* along the path of a photon. Finally, the ex- 
pression for the optical depth of the disk radiation fields to 
the TeV photons is written as 



T Tcv(^TeV,0Tev) - 



^77 (^t 



n v s ,n)-^-g\AQAv s Al, 
chv* 



(38) 

where dl = e v Y 1 A,9' is the proper length differential with 
d^ defined to be differential of the affine parameter 
along the TeV photons' trajectory (see Equation (22)), dVL = 
sinOsdOsdcpf, is the solid angle element in the LNRF, and /2 
is the cosine of angle between the TeV photons and the soft 
photons. The cross-section of the two colliding photons with 
energy e TcV = hv TcV /m e c 2 and e s = hvs/m e c 2 is given by Gould 
& Schreder (1967) 



3<7T n 2\ 



(3-v 4 )ln 



1+v 
L^v 



-2v(2-v 2 ) 



(39) 



for pair productions, where er T is the Thompson cross-section 
and v is the velocity of electrons and positrons at the center of 
momentum frame in units of c related to e TcV and e s through 



d-v 2 ) = 



(l-/x)e s e T 



(40) 



The cross section depends on the energies of the interacting 
photons and their colliding angle. If the two colliding pho- 
tons run in parallel (p, = 1), their interaction disappears. The 
lOTeV photons mostly interact with soft photons of 0.05eV 
for header-to-header collisions (p, = -1). As shown below, in 
specific calculations, the energy of soft colliding photons will 
be modified due to energy shift and bending of their trajecto- 
ries caused by the GR effects. 

It should be pointed out that the GR effects on r TeV occur 
through three factors: (1) changing the global structures of 
RIAFs, as a result, causing dependence of the radiation fields 
on the spin parameter; (2) modifying the observed spectrum 
from the RIAFs; and (3) bending the trajectories of TeV pho- 
tons. For different spins of black holes, the global structures 
of RIAFs become significantly different starting ^100 R g in- 
ward. Once considering the second influence, the observed 
spectrum is dependent on the viewing angles. The third influ- 
ence is nonnegligible unless the trajectories of TeV photons 
closely approach the black hole (~40 R g ). These three factors 
jointly influence the final optical depth to TeV photons. 

4. NUMERICAL RESULTS 

First, we numerically solve the GR RIAF equations for its 
structure spanning a large space of parameters. The results are 
consistent with the numerical solutions in Manmoto (2000). 
Second, we calculate the intrinsic spectrum for each case of 
RIAFs with different accretion rates and spins. Third, the ray 
tracing method is employed to get the spectrum in an ob- 
server's frame. Last, the optical depth to lOTeV photons is 
calculated along their path from (Rj e y, Qj e y) to the observer's 
sky at (oo,0 obs ). The constants of motion A Te v and Q TeV are 
obtained according to Equation (26), where the impact a is set 
zero corresponding to the shortest path and (3 is determined by 
solving the geodesic equations. The properties of structures 



> 




FlG. 3. — Spin dependence of r TcV for parameters m = 2.5 X 10 3 and 
©obs = ©TeV = 30° at different Rj e y. 

and emergent spectrum of RIAFs can be found in Manmoto 
(2000). Here we focus on the optical depth of TeV photons. 

Hereinafter we scale accretion rate in the units of the Ed- 
dington rate M = rhM^, where M^di = ^GM,m p /riaic and 
rj = 0.1 is the radiative efficiency. We set M. = 3.2 x 10 9 M Q 
throughout the paper, and ay = 0.1, /3d = 0.5, and 8 = 1.0 x 
10~ 3 in this section. 

4.1. spin dependence 

The black hole spins enhance the GR effects on the RIAF 
structures, the emergent spectrum, and the photon trajecto- 
ries. Detailed calculations show that both the surface density 
and the temperature of accretion flows at fixed radius increase 
with the spins, consequently, making the radiation more ef- 
ficient with the spins. In light of the drag of the frame, the 
trajectories of TeV photons will be elongated and the prob- 
abilities of pair productions are thus enhanced. These two 
factors together increase t TcV with the spins. However, the GR 
influence on the photon trajectories appears within a radius 
~10 R g and, therefore, is less important at relatively large ra- 
dius. 

Figure 3 shows the dependence of r TcV on the black hole 
spins for given parameters m = 2.5 x 10~ 3 and D b s = Qjev = 
30° at four emission radii of TeV photons. Obviously, r TeV 
increases with the spins, and its dependence on the spins 
is more tight at smaller radius because the GR effects be- 
come more significant. We can find that r TcV > 1 for a > 0.7 
at Rj e y = 15 R g . On the other hand, the TeV photons from 
^Tev = 6/? g cannot escape from the radiation fields of the GR 
RIAFs for all spin ranges since r TcV ^> 1 . 

4.2. Rj e \ dependence 

In Figure 4 we show r TcV as a strong function of /?iev, the 
radius of the location of TeV photons, by fixing m = 2.5 x 10" 3 
and 6 b s = Oiev = 30°. The radiation fields from RIAFs be- 
come more and more intensive when /?xev approaches the 
black hole, leading to a deep dependence of t TcV on /?iev- It 
is difficult to get an analytical formulation of the dependence, 
but Figure 4 shows the details of the dependence numerically. 
The radiation fields from the RIAFs have different spatial dis- 
tributions with the spins since the GR effects cause their struc- 
tures distinct. As explained in Section 4.1, the spins enhance 
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R t e y/ R g 

FIG. 4. — i?T e v dependence of T TeV for parameters a = 0.998, m = 2.5 X 10~ 3 , 
ande obs = e T eV = 30°. 

the radiation fields. If we fix the energy budget released from 
the RIAFs, the radiation fields should be more compact for the 
larger spins, causing a steeper T TeV — 7?TeV relation. The /?TeV 
dependence of t TcV can be applied to estimate the spins. 

4.3. bs and 0T e v dependence 

The dependence of t TcV on 9 b s and Oxev is caused by the 
anisotropy of radiation fields from the disk and the depen- 
dence of the cross-section on the colliding angle in particu- 
lar. Generally speaking, the number density of soft photons is 
mostly contributed from the inner region of the accretion disk. 
At larger Rj e y, the path of these photons that reach the inter- 
acting point with the Te V photons approximately has the same 
polar angle with 0j e v- Since when the two colliding photons 
run in parallel, their interaction disappears. We can see r TeV 
reaches its minimum at 9 b s = 0TeV, whereas at smaller Rj e y, 
the TeV photons are embedded in the radiation fields from the 
inner region of the accretion disk and interact with soft pho- 
tons from all directions and, consequently, it is difficult to give 
a general explanation to the dependence on b s and Oxev- 

We present the t TcV -dependence on 8 b s and Oj e v in Figures 
5 and 6, respectively, for the fixed parameters m = 6.3 x 10~ 3 
and a = 0.998 at different Rj e y. From Figure 5 we find that 
(1) for TeV = 0°, r TcV monotonously increases with 6 bs, 
and (2) there is a minimum t TcV for relatively large /? TeV at 
o bs = ©TeV- This is confirmed by the case of Oxev = 45°. 
Figure 6 shows that the t TcV -dependence on 9x e v has the sim- 
ilar properties as that on o b s . 

4.4. m dependence 

The t TcV -dependence on accretion rates is more straightfor- 
ward in light of the changes of number density of soft pho- 
tons and SED from the RIAFs. Figure 7 illustrates how r TcV 
changes with accretion rates for fixed parameters a = 0.998 
and <5> bs = 6xev = 30°. We find that t TcV is extremely sensi- 
tive to accretion rates as t TcV oc m 2 5-5 °, and the power index 
tends to be flat at larger /? TeV . This strongly indicates that TeV 
photons can be violently diluted by minor changes in accre- 
tion rates. In this sense, intensive 7-ray emission cannot be 
detected in variable sources. This has important implication 
in observations. 



This strong dependence can be explained in the following 
ways. Firstly, RIAFs have properties 

E ~ot, B~m l/1 , z/ syn <~ 7^m 1/2 , r es ~ m, (41) 

where B is the magnetic field, v sya is peak frequency of syn- 
chrotron emission, j e = kT e /m e c 2 is the Lorentz factor of hot 
thermal electrons, and r es is the optical depth for Thompson 
scattering. The synchrotron emission power approximates as 
<~ 7 2 £o« ~ 7e?w 2 and the bremsstrahlung emission power 
<~ m 2 . In our calculations, we find that T e oc rh q and hence 
7 e oc m q , where q <~ 1/3. To understand the t TcV -dependence 
on m, we have to know the energy of soft photons interact- 
ing with TeV photons. Figure 8 shows the contribution to 
t TcV for soft photons with different frequencies from an RIAF 
model with m = 4.0 x 10" 3 , a = 0.998, 9 obs = e Te v = 30°, and 
^Tev = 5/? g . We find that u s ~ 10 14 5 Hz, which corresponds 
to the first Comptonization peak or the second determined by 
the parameters of the RIAFs. Since RIAFs are optically thin, 
most of the synchrotron photons escape and a small fraction 
(<~ r es ) are Compton-scattered by the hot electrons, the en- 
ergy densities Uc, oc -fjm 3 oc m 3+2q and Uc 2 oc -f 2 m 4 oc rh A+2q 
for the first and second Comptonization, respectively. This 
determines the sensitivities of r TcV -dependence on m as 

f U c , cxm 3+2 «-m 3 ' 7 , 

T « K {t£ (42) 

For larger m, the bremsstrahlung emission will dominate the 
synchrotron emission and the Comptonization, even at the 
lower frequency. This leads to a flatter power index r TeV oc 
m 2 5 . These are nicely consistent with the numerical results as 
shown in Figure 7. 

We would like to point out that the very sensitive depen- 
dences of t TcV on accretion rate as t TcV oc m 2,5 " 5 evidently 
indicate r TcV >• 1 for the standard accretion disk, immedi- 
ately drawing a conclusion that the TeV photons cannot es- 
cape from the vicinity of SMBH fueled by the standard accre- 
tion disk. This is confirmed by the work of Zhang & Cheng 
(1997). 

5. APPLICATION TO M87 

TeV emission with a rapid variability (<~ 2 days) has been 
detected by the HESS in the giant elliptical galaxy M87, giv- 
ing clear evidence for a compact TeV emission region in the 
immediate vicinity of the black hole (Aharonian et al. 2006) 
and thus providing an opportunity to constrain the spin of the 
central SMBH. It is well known that M87 hosts an SMBH 
withM. = (3.2±0.9) x 10 9 M Q (Harms et al. 1994; Macchetto 
et al. 1997) at a distance of <~ 16 Mpc. Table 1 summaries the 
up-to-date observations from radio to X-ray band as shown in 
Figure 9 and 10. We exclude the data with low resolution so 
as to minimize the contamination from the knots in M87. The 
SED of M87 gives clear evidence for the RIAFs at work in its 
nucleus (Reynolds et al. 1996; Di Matteo et al. 2003). 

5.1. Accretion rates of the SMBH 

The radiation fields around the black hole can be quanti- 
fied from the spectral fitting of M87. Before spectral fitting, 
it is useful to understand to what extent the different parame- 
ters influence on the structures of the RIAFs and hence their 
radiation fields. Small «d indicates the low efficiency of an- 
gular momentum transfer, leading to the low radial velocity 
but high surface density; {3d represents the fraction of pres- 
sure contributed by the gas. As to lower /3<j, the structures 
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are changed through two ways: one bt enhancing the mag- 
netic field, and the other by increasing the temperature but 
reducing the surface density; m determines the overall peaks 
of the spectrum; the parameter 5, representing the fraction 
of the viscous dissipation that heats electrons, plays an im- 
portant role in determining the temperature of the electrons 
and, therefore, in the energy boost of photons after Compton- 
scattering, i.e., S mainly determines the frequency location of 
the Comptonization bumps. The value of S is still open to 
question at present since some nonthermal mechanisms (e.g., 
Kolmgorov-like turbulent cascades and collisionless shocks) 
are highly unclear, which may significantly change the value 



of (5 (Narayan & Yi 1995). We treat 5 as a free parameter in 
our calculations. In the spectral fit, we set ad and /3d a priori 
to reduce the freedom. We use the GR RIAF model to fit 
the multiwavelength spectrum of M87. We firstly fit the ob- 
servation data used in Paper I (hereafter Case I). We find that 
ad = 0.025 and (id = 0.2 give a good fit to the data, m and 
S are adjusted with different spins. We furthermore take into 
account the HST observations (Sparks et al. 1996; Maoz et al. 
2005) and redo the fit (hereafter Case II). We set the character- 
istic value of ad = 0.1 and (3d = 0.5 which means equipartition 
between the gas and magnetized fields. We present the spec- 
tral fit in the left panels of Figures 9 and 10 for Case I and 
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FIG . 7 . — m dependence of r TcV for parameters a = 0.998 and o bs = ©TeV = 
30° . From the left to the right, the lines correspond to the results for i?xeV = 
3,4,5,6,10, 15, 20, 30, 40J? g , respectively. We can find that t TcV ~ m 2 5 ' 5 . 
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FIG. 8. — Contribution to r TcV from different frequencies of soft photons 
from the RIAFs with parameters m = 4.0 X 10~ 3 , a = 0.998, obs = ©TeV = 
30° andR T eV =5R g . 

Case II, respectively. The fit parameters are listed in Table 3. 
The main differences between the two cases are the locations 
of the multi-Comptonization bumps. In terms of the energy 
dissipation to heat the electrons, the parameters ay and 6 can 
be absorbed into one parameter (see Equation (11)). We note 
that Joy of Case II is larger by factor ^4 to that of Case I, 
which leads to the difference by factor ~4 on the location of 
the first Comptonization bump. 

The accretion rates obtained from fitting of the SED in M87 
can be examined independently. Di Matteo et al. (2003) 
give the Bondi accretion rate of the nucleus of M87 using the 
Chandra observation as an upper limit of the accretion rate, 
M < M Bo ndi = 1 -6 x 10~ 3 M E dd- The accretion rate given by the 
RIAFs fitting in the current paper is consistent with this upper 
limit. 

We list published literature that gives an estimate of jet 
power in M87 in Table 2. In regard to the inevitable uncer- 



TABLE 1 

Summary of data for the nucleus of M87. 



Frequency 


vF v 


Resolution 


Ref. c 


Obs. 


u(Hz) 




















5.0 X 10 9 


0.1" 


0.7 


1 


VLBI 


2.2 X 10 10 


0.48 a 


0.15 


2 


VLBI 


1.0 x 10" 


8.7" 


0.1 


3 


VLBI 






IR 






Z.o X 1U 


AA A 1 "> ^ 
40.4 zt Z.J 


4oU 


1 


Gemini 






Optical-UV 






r A w 1A14 

6.0 X 10 


61.7 ± 12. j 


22 


6 


| ■ f \ f ' r7fT n 

r(JC, HS1 


7 n v in 14 


zu.u 


22 


A <Z 


ruo, no i 


8.1 x 10 14 


37.2 ±7.4 


22 


6 


FOC, HST 


9.1 x 10 14 


16.0 ±1.6 


28.4 


9 


ACS, HST 


1.2 x 10 15 


14.0 ±1.4 


28.4 


9 


ACS, HST 


1.3 x 10 15 


20.4 ±4.0 


22 


6 


FOC, HST 


1.9 x 10 15 


12.9 ±2.5 


22 


6 


FOC, HST 


2.0 x 1() 15 


16.2±3.2 


22 


6 


FOC, HST 


2.4 x 1() 15 


38.9±7.8 


22 


6 


FOC, HST 






X-ray 






2.4 X 1() 17 


8.0±0.2 b 


500 


8 


Chandra 



a No mention of error-bars. b The X-ray spectral index a x = -1 .23 ± 0.04, 
defined as F v oc . c References:- (1) Pauliny-Tofh et al. 1981; (2) 
Spencer & Junor 1986; (3) Baath et al. 1992; (4) Harms et al. 1994; (5) 
Reynolds et al. 1996; (6) Sparks et al. 1996; (7) Perlman et al. 2001; (8) 
Di Matteo et al. 2003; (9) Maoz et al. 2005. 



TABLE 2 

The jet power from the published 
literature 



L//erg s 1 


Ref. 


~ 10 44 


Bicknell & Begelman (1996) 


2 x 10 43 


Reynolds et al. (1996) 


~ 10 44 


Owen et al. (2000) 


3 x 10 42 


Young et al. (2002) 


~ 10 44 


Stawarz et al. (2006) 


5 x 10 43 


Bromberg & Levinson (2008) 



tainties of the estimates, the jet power in M87 is in a range of 
1 42 — 1 44 erg s" 1 . There are a number of theoretical calcu- 
lations of jet power driven by BZ (Blandford-Znajek) or BP 
(Blandford-Payne) mechanisms. For example, Meier (2001) 
developed a hybrid jet model by combining the two processes, 
in which the jet power is given by 

^Kerr-ADAF _j 

=0.1 (a^f) (0.14/ 2 + 0.74/j + ;V, (43) 

where / = fi / ^o,ny, j = cJ/GM 2 , g = B^o/B^ny, and L acc = 
O.lMc 2 (see Meier 2001 for details). Taking / = ; = g ~ 1 
and a = 0.1, the jet power will be Lj w 0.6L acc . Indeed, if 
considering the field-enhancing shear by frame-dragging ef- 
fects (Meier 2001, Nemmen et al. 2007), which are neglected 
in the self-similar solutions, we may have Lj > L acc , corre- 
sponding to extracting the rotational energy from the black 
hole (the Penrose mechanism). The accretion rates obtained 
in our paper m ~ 10" 4 are able to produce the jet power 
Lj > L acc = 3 x 10 43 erg s" 1 , which generally satisfies the jet 
energy budget from observations. 

5.2. Spin of the SMBH 
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FIG. 9. — Left panel: The spectral fit of M87 for Case I (see Table 2). The data from radio to X-ray band are summarized in Table 1 and the TeV data are from 
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TABLE 3 
The fit parameters (0 o b s = 30°). 





a 


0!d 




S 


n) 


R(r TeV = V>/R t 




0.0 


0.025 


0.20 


0.39 


2.0 X 10~ 4 


28 


Case I 


0.8 


0.025 


0.20 


0.18 


1.8 X 10~ 4 


20 




0.998 


0.025 


0.20 


0.09 


1.4 x 10~ 4 


15 




0.0 


0.1 


0.5 


0.35 


1.0 x 10"* 


28 


Case II 


0.8 


0.1 


0.5 


0.18 


9.0 x 10- 5 


20 




0.998 


0.1 


0.5 


0.11 


7.5 X 10~ 5 


16 



We calculate the optical depth to lOTeV photons and show 
the results in the right panel of Figures 9 and 10, in which we 
set o b s ~ 30° according to the VLB I observation of the jets 
in M87 (Bicknell & Begelman 1996) since the jets generally 
align at the axis of the accretion disk. For both cases, the 
resultant r TeV -^T e v relation becomes steeper with the spins. 
If we define the transparent radius R c as the radius at which 



r TcV = 1, we find R c 
in Case I, respectively, and R, 

II. The variability of TeV emission at a timescale of ^2days 
constrains the emission region Rj e y < 2QR g (Aharonian et al 



15R g and 28fl g for a = 0.998 and a = 
l6R g and 287? a in Case 



2006; Paper I). We find that a black hole with a spin of a TeV = 
0.8 for both cases leads to R c = 2QR g . To avoid being optically 
thick to lOTeV photons, it requires r TcV < 1 at R Te y = 20R g . In 
this sense, the spin of SMBH in M87 should be a > 0.8 as 
shown in Figures 9 and 10. 

We have to point out that the quality of the current data 
does not allow us to ascertain which fit is more practical to 
describe the accretion flows in M87. Future multiwavelength 
observations with high spatial resolution will provide stronger 
constraints on the spectrum and hence on the spin of the black 
hole in M87. 

6. DISCUSSIONS 

In this paper we fix the horizon of the black hole as the 
inner edge of the accretion flows. However, Krolik & Haw- 
ley (2002) show that the inner edge of the accretion flows is 
dependent on the accretion rate and is time- variable. To inves- 
tigate this effect on the optical depth, we compare the optical 
depth to lOTeV photons for the cases with inner edges set at 
the horizon and the last stable orbit, respectively. Figure 1 1 
shows the results with parameters a = (the horizon is 2R g 
and the last stable orbit is 6fl g ), a d = 0.1, /3 d = 0.5, 8 = 0.35, 
m = 1.0 x 10~ 4 , and Oxev = ©obs = 30°. We can clearly see 
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FIG. 1 1 . — Influence of the inner edge of the accretion disk on the spectra and r TcV . The parameters are a = 0.0, = 0.1, /3d = 0.5, <5 = 0.35, m = 1.0 X 10 -4 
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From Figure 1 1 that the inner iage lias little influence on the and 



optical depth of the TeV photons. We presume that the loca- 
tion of TeV origination should be close to the horizon of the 
black hole. This provides strong constraints on the radiation 
fields near the black hole horizon. With the detailed calcula- 
tions presented in this paper, there are still some aspects to be 
included in the future work. 

(1) The size of the TeV source. We simply assume a point- 
like TeV emission source and neglect its size. Once consider- 
ing spatial distribution of the source, r Ti . v should depend on the 
location of TeV photons since they have different trajectories. 
We here only focus on the minimum optical depth; however, 
the results in the present paper are viable for the current in- 
struments of TeV detection. Future work on the size effects 
could produce interesting results on variabilities of TeV pho- 
tons, including the profile of light curves and the time lag be- 
tween the different TeV photons. This may provide a mapping 
of the TeV source so that the radiation mechanism will be fi- 
nally discovered (Bottcher & Dermer 1995; Levinson 2000; 
Neronov & Aharonian 2007; Rieger & Aharonian 2008). 

(2) The motion of the TeV source. For TeV photons with 
specified trajectory, their optical depth is independent of the 
properties of the TeV source. However, the beaming effects 
caused by relativistic motion of the TeV source will help to 
avoid pair-production absorption of TeV photons, in regard to 
the beaming of the intensity of TeV photons along the direc- 
tion of the relativistic motion, and the blueshift of the energy 
of TeV photons by comparing with that in the source frame. 
This will modify the observed spectra of TeV emission. Inter- 
estingly, the present procedure modified by including beamed 
TeV photons from relativistic jets can be applied to blazars, 
which are being powered by the standard accretion disks (in 
flat spectrum radio quasars) and the RIAFs (in high-frequency 
peaked BL Lacs, such as Mkn 421 and Mkn 501) (Wang et al. 
2003). 

(3) The vertical structure of the RIAFs. An exponential pro- 
file of the density has been adopted in the vertical direction of 
the RIAFs, which means that most soft photons are emanat- 
ing from the midplane. Future three-dimensional simulations 
of the accretion flows will give a more realistic description of 
the RIAF structures. In this sense, the optical depth of the TeV 
photons should be calculated in a more sophisticated way. 

(4) The time-dependent trajectories of TeV photons. Fu- 
ture detectors with large area will receive more TeV photons 



ana be able to give more details of spectral variabilities. This 
needs a consideration on the time-dependent trajectories of 
TeV photons in light of r TeV (time-dependent pair productions 
when t TcV > 1), for example, the periodic rotation of the TeV 
source around the central SMBH. This will definitely provide 
much more information about the innermost region near the 
black hole's horizon. 

In other words, further detailed theoretical work is needed 
for stronger limits on SMBH spins from future observations. 
The current calculations based on the GR RIAF model pro- 
vide valuable constraints on the target fields of TeV photons 
and therefore on the spins. 

7. SUMMARY 

The optical depth to energetic TeV photons, which are im- 
mersed in the radiation fields from radiatively inefficient ac- 
cretion flows, has been calculated in detail by including all the 
GR effects. We investigate the dependence of optical depth 
on the spins (a), accretion rates (m), viewing angles (0 bsX 
and location of the TeV photons (Rj e y, Qj e y). We find that 
the optical depth is more sensitive to Rj e y than to 9 b s and 
©TeV- One of the most interesting results is that r TcV strongly 
depends on the accretion rates as r TeV oc to 2,5-50 . 

Applying the dependence of optical depth on Rj e \ to con- 
strain the spin parameter in M87, wherein the RIAFs are ex- 
pected to be at work, we find that the observed TeV photons 
detected by HESS can escape from the radiation fields from 
the RIAFs with spin a > 0.8. Future observations of the phase 
II HESS with threshold energy one order higher than that of 
the phase I may discover more M87-like objects (e.g., low- 
luminosity AGNs) with TeV emission. Hopefully, we then 
will have a sample for the statistic sample of spins of SMBHs. 
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APPENDIX 
PRELIMINARIES OF GR NOTATIONS 



We adopt geometrical units (G = c= 1) throughout the Appendix, where G is the gravitational constant and c is the light speed. 
We use the Kerr metric in Boyer-Lindquist coordinates (f , r, <fi, 6) 



with 



and 



ds 2 = -e 2u At 2 + e 2 ^ (d<f> - udt ) 2 + e 2 ^ 1 dr 2 + e 2 w d(9 2 



lv _ SA ^ _ sin 1 9A 

e A ' C E ' 

i, E t, „ 2M.ar 

e 2/« g 2 ^=E, uj= ^— , 



A = r 2 -2M.r+a 2 , £ = r 2 + a 2 cos 2 8, A = (r 2 + a 2 ) 2 -a 2 A sin 2 ( 



where a = J/M, is the specific angular momentum and M. is the black hole mass. The event horizon lies at 

r h =M. + (M.-a 2 ) 1 / 2 . 

The angular frequencies of the corotating (+) and counterrotating (-) Keplerian motions are 



fi± = ± 



M, 



1/2 



r 5 / 2 =F aM, 



1/2' 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 

(A6) 



To describe the motions of the accretion flows or photons in Kerr metric, we employ three reference frames (Gammie & 
Popham 1998). The first is the LNRF, an orthonormal tetrad basis carried by observers who live at constant r and 8, but at 
(f> = Lot + constant; the second is the CRF, whose coordinate angular velocity is ft, here ft is the angular velocity of the accretion 
flows. The last is the LRF of the accretion flows. The basis vectors for the LNRF are 



^ = ^(1,0,0,00, 
^ = ^(0,1,0,0), 

he)' 



e? m = e-^-(0, 0,1,0), 



e»=e^(0, 0,0,1), 



(A7) 
(A8) 
(A9) 
(A10) 
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and for the LRF are 

<) = (7r7*e"", 7r^"" 1 . 0, 7,70^0), 
<) = (7r7^"", 7r«"" 1 , 0, >7^ r <r^), 
e f e) = (0, o, 0), 

< 0) =[70^, 0, 0, 7 ^we-' + e-*)], 

where (3 r is the radial velocity of the accretion flows in the CRF with ~/ r = (\-0f)~ x l 2 , (3$ 
velocity of the CRF with respect to the LNRF with 7 = (1 - /3^)" 1/2 . 



(All) 
(A12) 
(A13) 
(A14) 

e^~ u (Q-u) is the physical azimuthal 



